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Abstract:--The paper presents the estimation of lateral aerodynamic parameters using neural based (Neural-Gauss-Newton) 
method from real flight data of Hansa-3 aircraft. The conventional methods require exact model postulation whereas the 
Neural-Gauss-Newton method is an algorithm that utilizes Feed Forward Neural Network and Gauss-Newton optimization 
to estimate the parameters and does not require a priori postulation of mathematical model or solution of equations of 
motion. The results obtained in terms of lateral-directional aerodynamic parameters were reasonably accurate to establish 
neural method with an additional advantage of non-requirement of a priori aerodynamic model. 
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I. INTRODUCTION 

The parameter estimation [1-10] from real flight data is a routine task for many aerospace organizations. Parameter 
estimation is the process of determining the best possible estimates of the parameters occurring in the model used to 
represent a system. Although reasonably accurate parameters can be obtained through analytical predictions and wind tunnel 
testing, the parameter estimation using flight data help to enhance the confidence in the estimates significantly. Designing 
optimal controls and autopilots, expansion of flight envelopes, updating simulators and verification of overall aircraft 
performance are some of the uses of parameter estimation. The conventional methods (Output Error method and its variants 
[2-5] such as Least Squares (LS) & Maximum Likelihood (ML) methods) used for identification assume the model to be 
exact. Most of the estimation methods find it difficult to handle flight data having reasonable amount of process noise. The 
problems of handling process noise and requirement of any mathematical model have been successfully dealt with in the 
present work by using an estimation method that utilizes Feed Forward Neural Network [6-10] and Gauss-Newton 
optimization that maps the input/output measurements of the system directly and has the capability to model any nonlinear 
continuous function without going into the physical details of the network. The method named as Neural Gauss-Newton [8- 
10] method was proposed by Peyada [8]. The neural based method [Neural Gauss-Newton (NGN)] has been shown to 
adequately estimate lateral-directional aerodynamic parameters from lateral-directional real flight data. The paper presents 
the generation of flight data, data-compatibility check, aerodynamic model, parameter estimation and concluding remarks. 

II. GENERATION OF FLIGHT DATA 

A flight test program using the Hansa-3 [Fig. 1] aircraft, an in-house fully instrumented research aircraft, was 
conducted at the Flight Laboratory, IIT Kanpur to gather the real flight data with the help of a data acquisition system. An 
onboard measurement system installed on test aircraft provided the measurements using dedicated sensors for a large 
number of signals such as aircraft motion variables, atmospheric conditions, control surface deflections etc. The 
measurements made in flight were recorded on board at a sampling rate of 50 Hz using a suitable interface with a standard 
Laptop computer. The three sets of lateral-directional flight data were acquired by executing the aileron/rudder control inputs 
during flight tests. The three lateral-directional flight data sets nomenclatured as HLD1, HLD2 & HLD3 (Where H and LD 
refer to Hansa-3 and Lateral-Directional respectively) are processed and presented graphically [Figs. 2-4] in terms of the 
motion variables. These figures present the variation of lateral-directional motion variables such as angle of sideslip (/?), roll 
angle (ji), yaw angle (y/), roll rate (p), yaw rate (r), linear acceleration (a T ) along y-axis and velocity (V) pertaining to doublet 
aileron and/or rudder (S a and/or S r ) control inputs. 




Fig. 1 The Hansa-3 research-aircraft 
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Fig. 2 Lateral-directional flight data: HLD1 
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Fig. 3 Lateral -directional flight data: HLD2 
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Fig. 4 Lateral-directional flight data: HLD3 



It can be seen that the trim angle for aileron and rudder are approximately zero and -7 degrees, respectively. The 
values of y/, p, r and a y corresponding to trim condition can be observed to be zero. The trim value of fl can be observed to be 
approximately zero. The velocity was kept around 32 ms"'. It can also be observed that the lateral variables (<f>, p) were 
affected when the ailerons were deflected from trim condition whereas the rudder deflection affects the directional variables 
(fi, r). A positive increase in aileron deflection from trim condition results in negative roll rate and reduction in bank angle 
whereas a positive increase in rudder deflection results in negative yaw rate and positive angle of sideslip 
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III. DATA COMPATIBILITY CHECK 

The data compatibility check, which is also known as flight path reconstruction [4], is an integral part of aircraft 
parameter estimation. The recorded real flight data are mainly corrupted by systematic errors like scale factors, zero shift 
biases and time shifts. These errors introduce data incompatibility; for example, data incompatibility would exist by way of 
the measured incidence angles not being in agreement with those reconstructed from the accelerometer and rate gyro 
measurements. The main aim of a data compatibility check is to ensure that the measurements used for subsequent 
aerodynamic model identification are consistent and error free (as far as possible). 
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Fig. 5 Data Compatibility Check: HLD1 



25 



30 



5 



15 



10 
Time, s 

Fig. 6 Data Compatibility Check: HLD2 
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The data compatibility check was carried out on both flight data sets pertaining to lateral-directional case using 
observation equations and the ML method. Figs. (5-7) present the measured and computed response of motion variables such 
as /?, (j>, a, 6 and V obtained during the data compatibility check from flight data pertaining to lateral-directional control 
inputs. It can be observed that the computed response compares well with the measured response for most of the motion 
variables (a, 8, P and <j>). A similar trend in the variation of velocity was also observed. The unknown parameter vector 
[equation (1)], representing scale factor and biases, was considered adequate for reconstructing lateral-directional dynamics 
of Hansa-3 aircraft. 



® = [Aa x Aa y Aa z Ap Aq Ar Kp A/?] 7 



(1) 



The Maximum Likelihood method was used to estimate the compatibility factors [equation (1)] from flight data 
pertaining to lateral-directional control inputs. Table 1 presents the values of scale factor (Kp) and biases 
(Aa x Aa y Aa z Ap Aq Ar Aft) estimated from lateral-directional flight data. The values mentioned in parentheses are the 
Cramer-Rao bounds suggesting the level of accuracy. It can be observed from Table 1 that the biases are negligible and scale 
factor is close to unity. Also, the values of Cramer-Rao bounds estimated along with compatibility factors are very low. The 
scale factor close to unity, negligible biases and very low values of Cramer-Rao bounds establish the high accuracy level of 
the data gathered during the flight testing. 

Table 1 : Data compatibility check: Lateral-directional flight data 



Factors — > 
Input J, 


Aa x 
(m/s 2 ) 


Aa y (m/s 2 ) 


Aa z (m/s 2 ) 


Ap 
(rad/s) 


Aq 
(rad/s) 


Ar 
(rad/s) 




Af} 
(rad) 


HLD1 


0.2478 


0.1827 


0.083 


-0.0006 


-0.0011 


0.0023 


1.0167 


-0.0122 




(0.0023) 


(0.0009) 


(0.0005) 


(0.0000) 


(0.0000) 


(0.0000) 


(0.0005) 


(0.0002) 


HLD2 


0.4557 


0.2303 


0.1257 


-0.0007 


-0.0008 


0.0018 


1.0804 


-0.0141 




(0.0027) 


(0.0015) 


(0.0006) 


(0.0000) 


(0.0000) 


(0.0000) 


(0.0102) 


(0.0003) 


HLD3 


0.1730 


0.1649 


0.0201 


-0.0002 


-0.0008 


0.0021 


0.9816 


-0.0041 




(0.0030) 


(0.0034) 


(0.0011) 


(0.0000) 


(0.0000) 


(0.0000) 


(0.0094) 


(0.0004) 


( ) Cramer-Rao Bounds 



IV. AERODYNAMIC MODEL 

The following lateral-directional state equations (simplified case) were used to postulate the aerodynamic model 
for the estimation of lateral-directional parameters. 



(3 = -r + | sin 0 
P = pV 2 S w c 
f = pV 2 S w c 
<j> = p 



pVSw 
2m 

QzCj+jxzCJ 
2(lxIz-4) 

QxCn+jxzCj) 

2(lxlz-liz) 



(2a) 
(2b) 
(2c) 
(2d) 



The side-force, rolling moment and yawing moment coefficient appearing in equation (2) are modeled as per equation (3). 



Cy = Cy a + Cy e P + C Yp g) + Cy r g) + Cy Sr S r 

Q = c lo + c h p + c lp g) + c lr g) + c ha s a + c hr s r 

C n = C no + C nf p + C np g) + C nr g) + C nSr 8 r 



(3a) 
(3b) 
(3c) 



The aim was to estimate the unknown parameter vector, © [equation (4)] using Regression methods from the 
lateral-directional flight data corresponding to the doublet aileron and/or rudder control inputs. 



® - [ C Y„ C Yp C Y p C Y r C Y Sr C l 0 Clp C l p C l r C l Sa C l Sr C„ 0 C nf) C„ p C„ r C ngr C n .j 



(4) 



V. PARAMETER ESTIMATION USING REGRESSION 

The Neural Gauss-Newton [8-10] method is an algorithm that utilizes the Feed Forward Neural Network (FFNN) 
and Gauss-Newton optimization to estimate the aerodynamic parameters. The neural model has been used to predict the time 
histories of motion variables at (fc+l) th instant given the measured motion variables corresponding to k th instant (where k = 1 
to n: n is the total number of discrete data points). For all the practical purposes of parameter estimation, this approach helps 
in building flight dynamic model (in restricted sense) using measured input-output data and does not require an a priori 
postulation of the mathematical model or solution of equations of motion. The algorithm of NGN method used to estimate 
the parameters with the help of block diagram given in Fig. 8 has been summarized below. 

a) As a first step, the measured flight data undergoes the data compatibly check. The measured motion variables are then 
transferred to the center of gravity for the further use during estimation process (Blocks 1 -3 of Fig. 8). 

b) The procedure followed for the neural network training using the FFNNs is explained in the blocks 3-8. 
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c) The block 9 checks the convergence criteria for the FFNNs training. Once the training is accomplished the trained neural 
model is used for the parameter estimation. 

d) Using the chosen aerodynamic model, the already trained neural model is used to calculate the system output Y(k). The 
input U(k) is constructed in the block 10 using aerodynamic model fed through block 15. The input U(k) is fed to block 5 to 
estimate the system output Y(k). 

e) The computed response Y(k) and dY(k)/d® (from block 11) are fed to the blocks 12 and 13 to update the aerodynamic 
parameter. The aerodynamic model is updated using new set of aerodynamic parameters in block 15. The computation 
through steps (d) to (e) is continued till the convergence criterion (block 14) is achieved. Once the convergence is achieved, 
aerodynamic parameters are estimated along with associated standard deviation. 
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Fig. 8 The schematic of NGN method 



The NGN algorithm was applied to the real lateral-directional flight data sets pertaining to doublet aileron -rudder 
control inputs to estimate the unknown parameter vector © [Eqn 4]. The residual error between the measured flight data 
Z(k + 1) and the estimated neural output Y(k + 1) was minimized to estimate the unknown parameter vector (©). The 
network parameters varied were the number of hidden layers (1-3), the number of neuron in the hidden layers (2-10), the 
learning rate (0.1-0.8), the momentum rate (0.1-0.8) and the number of iterations (100-4000). The network parameters finally 
chosen gave a good match between the true and the predicted values of the time histories of the variables. The final FFNN 
structure consisted of one hidden layer having the five neurons with a learning rate of 0.3 and the number of iterations equal 
to 2000. 

Figs. (9-11) present the measured and the trained response of the motion variables along y-body axis obtained during the 
process of training the neural model. It can be observed that the measured response of /?, </>, y/, p, r and a y matches well with 
the trained response. 

Figs. (12-14) present the measured and the estimated response of the motion variables such as angle of sideslip (/?), bank 
angle {</>), yaw angle (y/), roll rate (p), yaw rate (r) and acceleration (a y ) along y-body axis obtained during the process of 
parameter estimation from the real lateral-directional flight data sets using the NGN method. It can be observed that the 
measured response of /?, </>, <//, p, r and a y matches well with the estimated response for most of the motion variables for all 
the ten sets of flight data. 
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Fig. 9 Measured and Trained response: HLD1 
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Fig. 10 Measured and Trained response: HLD2 
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Fig. 11 Measured and Trained response: HLD3 



10 15 
Time, s 

Fig. 12 Measured and Estimated response: HLD1 



The parameter vector (©) given in Eqn 4 was estimated from the compatible real flight data by minimizing the cost 
function using NGN method. Table 2 presents the estimated lateral-directional aerodynamic parameters along with their 
Cramer-Rao bounds. The estimated parameters are compared to the wind tunnel estimates (col. 2). It can be observed that 
the estimated aerodynamic parameters such as C Yfs , C Yg , C lfs , C ls , C n/! andC ns are consistent and in close agreement 
with the wind tunnel estimates. The flight data sets gave consistent values of the estimated damping (C; and C„ r ) and the 
cross (C ir and C n ) derivatives (parameters). The obtained values of aerodynamic parameters such as C Y and C Yy were also 
consistent. However, the values of the estimated parameters such as C; 0 and C„ 0 are having opposite sign in contrast to the 
wind tunnel estimates but their value is quite small or negligible as desired for most of flight data sets. The aerodynamic 
parameter Cy 0 could not be estimated correctly. The estimated value of the parameter C Ys was consistent but on a higher side 
for most of the flight data sets. The value of the parameter C; also could not be estimated correctly. 
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Table 2: Parameter Estimation using NGN Method 



Parameters 


WT Value 


HLD1 


HLD2 


HLD2 




-0.013 


0.0168 


0.0298 


0.0154 




(0.0003) 


(0.0005) 


(0.0006) 




-0.531 


-0.8096 


-0.8975 


-0.6154 






(0.0063) 


(0.0076) 


(0.0156) 






0.1825 


0.0375 


-0.0004 






(0.0219) 


(0.0248) 


(0.0001) 






0.6598 


1.3262 


0.0002 






(0.0273) 


(0.0350) 


(0.0001) 




0.150 


0.2254 


0.3733 


0.2021 




(0.0031) 


(0.0044) 


(0.0049) 




0.0015 


0.0007 


0.0012 


0.0008 






(0.0001) 


(0.0002) 


(0.0002) 




-0.031 


-0.0290 


-0.0253 


-0.0394 






(0.0021) 


(0.0032) 


(0.0047) 


c t 

L v 




-0.2865 


-0.3056 


0.0000 






(0.0088) 


(0.0137) 


(0.0000) 






0.1610 


0.1815 


0.0000 




(0.0091) 


(0.0141) 


(0.0000) 




-0.153 


-0.1549 


-0.1760 


-0.1258 






(0.0029) 


(0.0046) 


(0.0054) 


^ 


0.005 


0.0202 


0.0252 


0.0218 






(0.0011) 


(0.0017) 


(0.0021) 




0.001 


-0.0076 


-0.0109 


-0.0042 




(0.0001) 


(0.0001) 


(0.0002) 


Cn, 


0.061 


0.0323 


0.0462 


0.0495 






(0.0021) 


(0.0029) 


(0.0056) 


Cn p 




-0.1548 


-0.1427 


0.0001 






(0.0072) 


(0.0096) 


(0.0000) 


Cn r 




-0.1431 


-0.1654 


-0.0001 






(0.0092) 


(0.0131) 


(0.0000) 


C n Sr 


-0.05 


-0.0727 


-0.0997 


-0.0442 






(0.0010) 


(0.0016) 


(0.0018) 
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VI. CONCLUSION 

The NGN method was used to model the lateral-directional aerodynamics using flight data of Hansa-3 aircraft. The 
following points were observed during the aerodynamic modeling in time domain. 

1 . The correct calibration is necessary before the flight test program to acquire flight data of good quality. 

2. Any estimation method (output error methods) will give reasonably good estimates of parameters, if real flight 
data is of good quality. 

3. The estimated compatibility factors established the high quality of flight data gathered during flight testing. 

4. The NGN methods resulted in satisfactory estimation of lateral-directional aerodynamic parameters. 

5. The reason for the deviation of some parameters being the insufficient information content in the data generated. 

6. Weak parameters are difficult to estimate. Therefore, this could be the second reason for non-estimation few 
parameters (weak parameters). 

The NGN method does not require any a priori postulation of mathematical model or solution of equations of motion and act 
as functional approximator which can collectively model any nonlinear relationship between the inputs and the outputs and 
thereby provide overall characterization of a system. 
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